Search results for "Hessian matrix"
showing 10 items of 25 documents
PDF reweighting in the Hessian matrix approach
2014
We introduce the Hessian reweighting of parton distribution functions (PDFs). Similarly to the better-known Bayesian methods, its purpose is to address the compatibility of new data and the quantitative modifications they induce within an existing set of PDFs. By construction, the method discussed here applies to the PDF fits that carried out a Hessian error analysis using a non-zero tolerance $\Delta\chi^2$. The principle is validated by considering a simple, transparent example. We are also able to establish an agreement with the Bayesian technique provided that the tolerance criterion is appropriately accounted for and that a purely exponential Bayesian likelihood is assumed. As a practi…
Impact of dijet and D-meson data from 5.02 TeV p+Pb collisions on nuclear PDFs
2020
We discuss the new constraints on gluon parton distribution function (PDF) in lead nucleus, derivable with the Hessian PDF reweighting method from the 5.02 TeV p+Pb measurements of dijet (CMS) and $D^0$-meson (LHCb) nuclear modification ratios. The impact is found to be significant, placing stringent constraints in the mid- and previously unconstrained small-$x$ regions. The CMS dijet data confirm the existence of gluon anti-shadowing and the onset of small-$x$ shadowing, as well as reduce the gluon PDF uncertainties in the larger-$x$ region. The gluon constraints from the LHCb $D^0$ data, reaching down to $x \sim 10^{-5}$ and derived in a NLO perturbative QCD approach, provide a remarkable…
Can we fit nuclear PDFs with the high-x CLAS data?
2020
AbstractNuclear parton distribution functions (nuclear PDFs) are non-perturbative objects that encode the partonic behaviour of bound nucleons. To avoid potential higher-twist contributions, the data probing the high-x end of nuclear PDFs are sometimes left out from the global extractions despite their potential to constrain the fit parameters. In the present work we focus on the kinematic corner covered by the new high-x data measured by the CLAS/JLab collaboration. By using the Hessian re-weighting technique, we are able to quantitatively test the compatibility of these data with globally analyzed nuclear PDFs and explore the expected impact on the valence-quark distributions at high x. W…
Interactive simulation of one-dimensional flexible parts
2006
Computer simulations play an ever growing role for the development of automotive products. Assembly simulation, as well as many other processes, are used systematically even before the first physical prototype of a vehicle is built in order to check whether particular components can be assembled easily or whether another part is in the way. Usually, this kind of simulation is limited to rigid bodies. However, a vehicle contains a multitude of flexible parts of various types: cables, hoses, carpets, seat surfaces, insulations, weatherstrips... Since most of the problems using these simulations concern one-dimensional components and since an intuitive tool for cable routing is still needed, w…
A Sequential Quadratic Programming Method for Volatility Estimation in Option Pricing
2006
Our goal is to identify the volatility function in Dupire's equation from given option prices. Following an optimal control approach in a Lagrangian framework, we propose a globalized sequential quadratic programming (SQP) algorithm with a modified Hessian - to ensure that every SQP step is a descent direction - and implement a line search strategy. In each level of the SQP method a linear-quadratic optimal control problem with box constraints is solved by a primal-dual active set strategy. This guarantees L^1 constraints for the volatility, in particular assuring its positivity. The proposed algorithm is founded on a thorough first- and second-order optimality analysis. We prove the existe…
Impact of CMS 5.02 TeV dijet measurements on gluon PDFs - a preliminary view
2018
We discuss the implications of the preliminary CMS dijet data from 5.02 TeV pp and pPb collisions for gluon PDFs of the proton and nuclei. The preliminary pp data show a discrepancy with NLO predictions using for example the CT14 PDFs. We find that this difference cannot be accommodated within the associated scale uncertainties and debate the possible changes needed in the gluon PDF. A similar discrepancy is found between the CMS pPb data and NLO predictions e.g. with the EPPS16 nuclear modifications imposed on the CT14 proton PDFs. When a nuclear modification ratio of the pp and pPb data is constructed, the uncertainties in the scale choices and in proton PDFs effectively cancel and a good…
Real symplectic formulation of local special geometry
2006
We consider a formulation of local special geometry in terms of Darboux special coordinates $P^I=(p^i,q_i)$, $I=1,...,2n$. A general formula for the metric is obtained which is manifestly $\mathbf{Sp}(2n,\mathbb{R})$ covariant. Unlike the rigid case the metric is not given by the Hessian of the real function $S(P)$ which is the Legendre transform of the imaginary part of the holomorphic prepotential. Rather it is given by an expression that contains $S$, its Hessian and the conjugate momenta $S_I=\frac{\partial S}{\partial P^I}$. Only in the one-dimensional case ($n=1$) is the real (two-dimensional) metric proportional to the Hessian with an appropriate conformal factor.
Hessian PDF reweighting meets the Bayesian methods
2014
We discuss the Hessian PDF reweighting - a technique intended to estimate the effects that new measurements have on a set of PDFs. The method stems straightforwardly from considering new data in a usual $\chi^2$-fit and it naturally incorporates also non-zero values for the tolerance, $\Delta\chi^2>1$. In comparison to the contemporary Bayesian reweighting techniques, there is no need to generate large ensembles of PDF Monte-Carlo replicas, and the observables need to be evaluated only with the central and the error sets of the original PDFs. In spite of the apparently rather different methodologies, we find that the Hessian and the Bayesian techniques are actually equivalent if the $\Delta…
Nuclear parton distribution functions with uncertainties in a general mass variable flavor number scheme
2020
In this article we obtain a new set of nuclear parton distribution functions (nuclear PDFs) at next-to-leading order and next-to-next-to-leading order accuracy in perturbative QCD. The common nuclear deep-inelastic scattering (DIS) data analyzed in our study are complemented by the available charged-current neutrino DIS data with nuclear targets and data from Drell-Yan cross-section measurements for several nuclear targets. In addition, the most recent DIS data from the Jefferson Lab CLAS and Hall C experiments are also added to our data sample. For these specific datasets, we consider the impact of target mass corrections and higher twist effects which are expected to be important in the r…
Analytic evaluation of the dipole Hessian matrix in coupled-cluster theory
2013
The general theory required for the calculation of analytic third energy derivatives at the coupled-cluster level of theory is presented and connected to preceding special formulations for hyperpolarizabilities and polarizability gradients. Based on our theory, we have implemented a scheme for calculating the dipole Hessian matrix in a fully analytical manner within the coupled-cluster singles and doubles approximation. The dipole Hessian matrix is the second geometrical derivative of the dipole moment and thus a third derivative of the energy. It plays a crucial role in IR spectroscopy when taking into account anharmonic effects and is also essential for computing vibrational corrections t…